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ABSTRACT 


The  well-known  deflection  coefficient  for  a  stationary  Gaussian  signal  in  ad¬ 
ditive,  independent,  stationary  Gaussian  noise  is  derived.  The  optimal  pre¬ 
detection,  or  Eckart,  filter  is  also  derived.  It  is  shown  that  the  least  detectable 
Gaussian  signal  of  specified  power  is  proportional  to  the  normalized  square  of 
the  noise  spectrum.  It  is  also  shown  that  maximally  detectable  Gaussian  signals 
of  specified  power  must  have  at  least  one  spectral  component  of  finite  power  and 
very  narrow  bandwidth. 
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X.  INTRODUCTION 


The  processing  stream  is  depicted  in  Fig.  1.  Pre-detection  and  post-detection 
filters  are  assumed  linear  and  time  invariant.  The  input  v(t)  is  assumed  real 
valued,  zero  mean,  Gaussian,  and  stationary,  so  x(t)  is  also  real  valued,  zero 
mean,  Gaussian,  and  stationary.  Due  to  the  nonlinear  square  law  device,  y(t) 
and  z{t)  are  real  valued  and  stationary,  but  neither  zero  mean  nor  Gaussian.  (If 
the  input  process  is  complex  valued  instead  of  real  valued,  as  assumed  here,  the 
Eckart  filter  and  corresponding  deflection  coefficient  differ  from  the  expressions 
given  below  by  a  constant.) 

Assuming  small  signal  to  noise  ratio,  the  variance  under  the  null  hypothesis 

{Ho  :  noise  only  is  present} 

is  approximately  equal  to  the  variance  under  the  alternative  hypothesis 

{Hi  :  signal  plus  noise  is  present}. 

The  difference  of  the  means  of  the  output  z(t )  under  these  hypotheses  is  divided 
by  the  standard  deviation  under  the  null  hypothesis  to  obtain  the  so-called 
deflection  coefficient,  denoted  by  d.  The  receiver  operator  characteristic  (ROC) 
curves  of  the  detector  are  completely  parameterized  by  d  if  Gaussian  signal  and 
noise  statistical  assumptions  hold. 

2.  DETECTOR  VARIANCE  UNDER  NOISE-ONLY  HYPOTHESIS 

The  variance  of  the  output  z(t)  is  computed  first.  Let  hi  (r)  denote  the  impulse 
response  function  of  the  post-detection  filter  normalized  so  that 

[  hi(r)  dr  =  1.  (2.1) 

J—OO 

The  output  z(t)  is  the  convolution  of  the  input  y{t)  with  hi  (r),  so  that 

z(t)  =  [  y(t -r)hi(r)  dr.  (2.2) 

J—OO 

Substituting  (2.2)  into  the  definition  of  the  autocorrelation  function  of  z(t)  gives 
Rz(r)  =  E[z(t)z(t  +  r)] 
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=  e\[  y(t  -  a)hi(a)  da  [  y{t  +  t  - /3)hi(j3)  d/3 

[y  — oo  J —oo 

/oo  roc 

/  £[y(f  -  a)y(f  +  r  -  /3)]hi(a)/ii(/3)  da  d/3 

-00  */  — OO 

/OO  TOO 

/  i^(r  -  /3  +  a)hi(a)hi(J3)  da  d/3,  (2.3) 

•oo  J  —oo 

where  the  autocorrelation  function  of  y(t)  is 

i?y(r)  =  £[y(t)y(t  +  r)] 

=  E[a:2(t)x2(t +  r)] 

=  ^(0)+2i?2(r),  (2.4) 

and  where  Rx(r)  =  E[x(t)x(t  +  r)]  is  the  autocorrelation  function  of  x(t).  The 
result  (2.4)  is  found  in  many  places;  for  example,  Thomas1  derives  it  from  the 
fourth  moment 

E[x  1X2X3X4]  =  E[XiX2]E[x3X4]  +  E[xiX3]E[x2Xi] 

+E[xiX4]E[x2X3]  —  2  X1X2X3X4 

for  general  quadrivariate  Gaussian  random  variables  by  setting  x\  =  X2  =  x(t) 
and  x3  =  x4  =  x(t  +  r)  and  using  the  fact  that  x(t)  is  zero  mean.  Since 
E[z2]  =  Rz{ 0)  and  E[z\  =  R^( 0),  the  variance  of  z  is 

Var[z ]  =  E[z2]-E2[z] 

=  ft(0)-J^(0) 

/oo  fOO 

/  Ry(—f3  +  a)hi(a)hi(j3)  da  d(3  —  Rl(0) 

-00  J — OO 

=  jT  {^(0)  +  +  a)}  hi(a)hi{J3)  da  d/3  -  Jg( 0) 

roo  roo 

=  2/  /  i?£(-/3  +  a)hi(a)hi(/3)  da  d/3,  (2.5) 

where  the  normalization  (2.1)  is  used  in  the  last  step.  The  result  (2.5)  is  easily 
evaluated  for  any  specified  post-detection  filter. 

The  simplest  low-pass  post-detection  filter  satisfying  the  normalization  (2.1) 
is  h\(t)  =  1/T  for  |t|  <  T/2  and  hi(t)  =  0  for  \t\  >  T/2.  It  is  assumed  that  the 
averaging  time  T  is  large  compared  to  the  correlation  time  of  x(t).  Substituting 
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into  (2.5  )  gives 


Var[z]  =  f  j_T&x(T)\l-Yj  dT 
2 

“  ?/.*«* 

=  (2.6) 

where  the  power  spectrum  Sx(u>)  of  x(t)  is  the  Fourier  transform 

Sx(u)  =  r  i?I(r)e-iWT  dr.  (2.7) 

J  — oo 

The  last  equation  in  (2.6)  is  Parseval’s  theorem  for  (2.7),  namely, 

/  i£(r)  dr  =  -^-J  S%(v)  du. 

J  — oo  Z7T  J  —oo 

Denoting  the  impulse  response  and  transfer  functions  of  the  pre-detection 
filter  by  ho(r)  and  H0(u>),  respectively,  gives 

5x(w)  =  |^0(a;)|251,(a;).  (2.8) 

Substituting  (2.8)  into  (2.6)  gives 

Var[z]  =  (2.9) 

The  standard  deviation  of  z(t)  is  the  square  root  of  (2.9). 

3.  DEFLECTION  COEFFICIENT 

The  difference  in  the  means  of  z(t)  under  hypotheses  Ho  and  Hi  is  now  com¬ 
puted.  Denote  the  correlation  functions  and  power  spectra  of  the  input  v(t)  un¬ 
der  hypotheses  Ho  and  Hi  by  Rg(r)  and  S*(u),  and  by  Rs+N  (r)  and  S^+N  (u), 
respectively.  Signal  and  noise  axe  statistically  independent,  so  their  power  spec¬ 
tra  add;  that  is, 

Ssv+N(u)  =  Ssv(u)  +  S-(u),  (3.1) 
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where  S*  ( u )  denotes  the  spectrum  of  the  signal  when  signal  only  is  present  in 
v(t).  The  pre-detection  filter  is  linear  and  time  invariant,  so  the  difference  in 
the  mean  of  the  output  z(t)  under  hypotheses  Ho  and  Hi  is 


A n  =  E[z\Hi]  -  E[z\Hq\ 

=  -«?(<>) 

=  i  r  ^  i  r  & 

Z7T  J—  oo  Z7T  J  — oo 

1  r°° 

=  (3.2) 


where  in  the  last  step  (3.1)  has  been  used. 

The  deflection  coefficient  for  Gaussian  random  variables  with  identical  vari¬ 
ances  is  defined  by 


d  = 


A  fj, 


(3.3) 


where  a  denotes  the  common  standard  deviation.  Because  of  the  low  signal  to 
noise  ratio  (SNR)  assumption,  the  standard  deviations  of  z(t)  under  Ho  and  Hi 
are  approximately  equal.  Substituting  (3.2)  and  (2.9)  into  (3.3)  gives 


j  lv^rjwi2s>)^ 


(3.4) 


If  the  signal  spectrum  is  unknown,  the  pre-detection  filter  is  sometimes  omitted, 
so  that 

(3.5) 

H5 

which  is  the  expression  used  for  the  conventional  processor. 


4.  ECKART  FILTER 


The  deflection  (3.4)  is  easily  maximized  in  any  specified  frequency  band,  say  T', 
for  which  the  noise  spectrum  is  bounded  away  from  zero;  that  is,  Sg  (u)  >  e  >  0 
for  u>  €  Tr.  For  u  €  write 


[^(a.)!2^^)  =  [|i?0(a.)|2^(^)] 


'SM‘ 
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Prom  the  Cauchy-Schwaxtz  inequality, 


[  |tf0(w)|2Sj(w)  du  < 


Substituting  (4.1)  into  (3.4)  gives  the  inequahty 


1 

2 


du 


(4.1) 


d  < 


1 


f 

[s»l 

2 

du 

(4.2) 


The  inequality  (4.2)  holds  for  all  H0(u).  The  upper  bound  is  independent  of 
H(j(cj),  so  equality  holds  in  (4.2)  if  and  only  if  equality  holds  in  (4.1).  The 
inequality  (4.1)  is  an  equality  if  and  only  if,  for  some  (real)  proportionality 
constant  k, 

\H0(u)\2  St  («,)  =  K  (4.3) 

for  all  u  €  The  deflection  d  is  invariant  to  the  scale  of  Ho(u),  so  k  is  set  equal 
to  one.  Solving  (4.3)  for  \H0(u)\2 ,  and  relabeling  it  H Eckart (w)j  gives 

(4.4) 

[^(w)] 

The  optimum  filter  (4.4)  is  often  called  the  Eckart  filter,  after  its  discoverer.2 
Substituting  the  optimum  filter  (4.4)  into  (3.4)  gives 


dEckart  — 


1 

-  I 

2 

du 

2 

7T  J' 

(4.5) 


The  last  result  shows  that  {dEckart)2  is  directly  proportional  to  the  averaging 
time  T  of  the  post-detection  filter  and  the  integral  of  the  square  of  the  SNR. 

Scholium.  Eckart ’s  original  derivation2  does  not  use  the  Cauchy-Schwartz 
inequahty.  Instead,  he  writes  a  candidate  expression  for  the  optimum  filter,  and 
then  proves  that  this  expression  is  indeed  optimal.  He  gives  no  indication  of 
how  he  initially  obtained  the  correct  filter. 
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5.  EXTREMAL  SIGNALS  OF  SPECIFIED  POWER 


5.1.  Least  Detectable  Signal 

Given  the  noise  spectrum  Sg  (w)  and  a  specified  signal  power,  say  Sp^er,  the 
least  detectable  signal  Sopt(u)  is  the  signal  for  which  the  deflection  coefficient 
dEckart  is  a  minimum.  Let  a  >  0  denote  the  SNR.  Then  S^t{uj)  is  the  solution 
of  the  following  calculus  of  variations  problem: 


subject  to  the  SNR  constraint 


Jpower 


S(u) 

Sff(u) 


U  S(w)  dw 


U  Sp  (w)  dw  S»  (a ;)  dw 

Let  S( w)  =  f(w)  Sg  (u),  where  f(w)  >  0,  and  let 

g(w)  =  r^rhr  -  °-  ( 

USp(v)  dw 

The  problem  (5.1)— (5.2)  written  in  terms  of  f(w)  and  G(w)  takes  the  form 


subject  to 


min  /  f2(w)  dw 

/(«)  h 

j  f(w)G(w)  dw  =  a. 


It  is  shown  below  that  the  solution  to  (5.4)-(5.5)  is  non-negative,  so  the  con¬ 
straint  f(w)  >  0  is  ignored  here  without  loss  of  generality.  The  appropriate 
Lagrangian  function  for  (5.4)-(5.5)  is 

£(/)  =  \J^f(u)du,-  A  (jf  dw  -  a)  .  (5.6) 

Let  6f(w)  denote  a  fixed  variation  of  f(w),  and  let  e  be  any  real  number.  Then 
the  variation  of  (5.6)  is 

A/  +  e8f)  =\/<Sf  [/(w)  +  ^/(w)]2  ^  ~  A  (_£  +  G (w)  dw  -  a)  . 

(5.7) 
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The  derivative  of  C(f  +  e6f)  with  respect  to  e  must  be  zero  at  e  =  0;  hence, 


dC(f  +  e6f) 


de 


«=o 


=  [f(cu)  -  X G(u>)]  8f{uj)  du  =  0,  (5.8) 


for  all  variations  8f(u).  For  (5.8)  to  hold  for  all  <5/(a>),  it  must  be  the  case  that 
f(oj)  —  XG(uj)  =  0.  Substituting  f(u>)  =  XG(u)  into  the  constraint  (5.5)  and 
integrating  over  u>  gives 


a  =  J  [AG(u;)]  G(u)  cLi  =  X  J  G?(u>)  dm. 


Solving  the  last  equation  for  A  and  substituting  into  f(u)  =  XG(uj)  gives  the 
optimal  solution 


fopt(w)  — 


aG(  u) 
U&iv)  dLJ 


>0. 


(5.9) 


Because 


&£{f  +  eSf) 


de 2 


e=0 


dhJ  >  0 


for  all  nonzero  variations  the  sufficient  conditions  for  the  solution  (5.9)  to 

be  a  minimum  hold;  thus,  fopt(w)  is  the  unique  global  minimum  of  (5.4)-(5.5). 
In  the  original  problem  (5.1)— (5.2) ,  the  signal  that  minimizes  the  deflection 
coefficient  dsckart  is 


Sopt{u)  — 


rfirwr-  w  mi2 

U  isv  M] 

^  /*[^(u;)]2  d/ 


(5.10) 


In  words,  the  spectrum  of  the  least  detectable  signal  is  directly  proportional  to 
the  normalized  square  of  the  noise  spectrum,  and  the  proportionality  constant 
is  the  specified  signal  power. 


5.2.  Maximally  Detectable  Signals 

If  noise  power  is  finite  and  its  spectrum  is  continuous  and  bounded  on  the  band 
\Er,  then  the  deflection  coefficient  d-Eckart  can  be  made  arbitrarily  large  for  any 
specified  signal  power  SpoWer  >  0.  For  example,  if  the  signal  spectrum  is 

rf  /  \  f  ^  Spower  i  <1  ^  U  ^  CL  6, 

c  U;  |  0,  otherwise, 
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where  a  and  e  >  0  are  such  that  the  interval  [a,  a  4-  e]  is  a  subset  of  then 

f  Se{u)  du  =  fa+e  S€(uj)  dw  =  Sp 

Jv  Ja 


J  power 


and 


L 


'SAu) 

Sf'fc') 


S 2  ra+£ 


ra+e 

/  [W] 

Ja 


-2 


du. 


(5.11) 


By  the  mean  value  theorem  of  calculus,  there  exists  a  point  ue  <E  [a,  a  +  e]  such 
that 

/°+t  [^(a;)]-2  do;  =  e [S? (««)]' “2  .  (5.12) 

Ja 

Substituting  (5.12)  into  (5.11)  gives 

1  2 


[ 

r  se(w)  l 

Jy 

S£(w). 

1  9s  1 

du  =  -  ~~^S!SL-  >  ± 


c  [5^(w«)] 


2  — 


Jpovjer 


max  Sf?  (w) 

Uef  v  v  ' 


(5.13) 


By  assumption,  the  noise  spectrum  is  bounded  on  the  band  'P,  so  the  maximum 
in  (5.13)  is  finite.  Thus,  dsckart  is  made  arbitrarily  large  by  choosing  e  sufficiently 
small.  In  words,  the  most  detectable  signals  are  those  whose  spectra  contain  at 
least  one  component  of  finite  power  and  very  narrow  bandwidth. 


6.  SUMMARY 

The  deflection  coefficient  d  for  a  Gaussian  signal  in  additive,  independent  Gaussian 
noise  is  derived.  The  optimal  pre-detection,  or  Eckart,  filter  is  also  derived.  The 
signal  of  specified  power  that  is  least  detectable  by  the  optimal  processor  is 
directly  proportional  to  the  normalized  square  of  the  noise  spectrum,  and  it  is 
obtained  using  a  calculus  of  variations  argument.  It  is  also  shown  that  the  most 
detectable  signals  of  specified  power  have  at  least  one  spectral  component  of 
finite  power  and  very  narrow  bandwidth. 
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Figure  1.  Processing  Stream 
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